Bose-Einstein Condensation of Molecular Hydrogen in Nanotube Bundles 
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We evaluate the effects of heterogeneity on the density of states of H2 molecules inside interstitial 
channels within bundles of carbon nanotubes. As temperature (T) falls, the density increases within 
those tubes having the greatest binding energy. At T ~ 10 mK, the molecules undergo Bose-Einstein 
condensation, exhibiting a singular heat capacity. 



The subject of adsorption within bundles of carbon 
nanotubes has received considerable attention recently, 
owing to both its fundamental interest and potential ap- 
plications (e.g. gas storage and isotope separation). One 
focus of the research is one-dimensional (ID) and quasi- 
1D phases of matter. These include condensing and crys- 
tallizing phases of buckyballs within tubes, He and H 2 
within the interstitial channels (IC's) between tubes and 
various gases withingrooves on the outside of the nan- 
otube bundles IliroilllHSHInlGIl Most anal- 
yses (including those of our group) have assumed that 
the tubes are identical and parallel, forming an ordered 
lattice. Real nanotube bundles, in contrast, consist of a 
disordered array of tubes with a distribution of radii. A 
logical question arises: how reliable are predictions that 
ignore such a variable environment? Shi and Johnson 
have recently shown that predictions incorporating het- 
erogeneity agree better with adsorption data than do the 
idealized models |l2^ . Stimulated by their work, we have 
explored the behavior of quantum fluids (He and H2) in 
such a nonuniform environment at low temperature (T). 
In this paper, we describe an intriguing result: Bose- 
Einstein condensation (BEC) of H2 molecules occurs as 
a consequence of the heterogeneity. This paper makes 
predictions about this phenomenon that are testable ex- 
perimentally. Similar behavior is expected for 4 He atoms. 
We note that a superfluid phase of para-H2 in confined 
geometries has been proposed, but experimental evidence 
of it is lacking thus far 13]. 

The occurrence and properties of the BEC transition 
are determined by the density of states N(E) of the con- 
stituent particles. Thus, an initial focus is the determina- 
tion of that function, derived from the single particle en- 
ergy spectrum. We assume that particles do not interact 
(based on the assumption of low density and weak inter- 
IC interactions) . A typical IC presents a highly confining 
geometry for the molecules if the bundle is close-packed. 
Because of this confinement (in the x-y plane, perpendic- 
ular to the tubes), the molecules have a large zero-point 
energy of motion. Theoretical values of this energy 
are of the order of 500 K. This result is consistent with an 
enormous isotopic heat difference at low coverage (which 



equals the binding energy difference) between D? and H2 
in a nanotube bundle observed experimentally |15| . That 
difference was found to be about 250 K, a factor w four 
larg er than the difference found on the graphite surface 
|l61| (due to less localization in that case). 

Our analysis assumes that molecules can move within 
the IC's in order to achieve chemical equilibrium with 
a coexisting vapor phase. At low T, because of the 
confinement, only the lowest energy state of transverse 
motion is excited. We call this state's energy i£ t (R), 
where R = (Ri, R2, R3) is a vector whose components 
are the radii of the tubes surrounding a particular IC. 
We have evaluated E t (TL), with interesting results. The 
calculation assumes that motion parallel to the IC is 
that of free molecules, so the total energy of a par- 
ticle with (z component of) momentum p is given by 
E(p,K) = E t (R) +p 2 /{2m). If the potential is corru- 
gated, one must replace the particle's mass m with its 
band mass 0, . The density of states for H2 is ob- 
tained by summing over the IC's present in the given 
sample and integrating over p: 

N(E) = S[E-E(p,K)} 

= ii (™) 1/2 [ E dE t -0L (i) 

Kit \ 2 ) J ^E^W t 1 ' 

Here g{E t ) is the density of states for the transverse 
oscillation problem and L is the length of the tubes. If 
the length is not constant, the variation can be included 
in g(E t ). 

The energy E t (R) is evaluated from the potential en- 
ergy F(r,R) of the molecule at position r. To com- 
pute V, we add contributions from the three neighboring 
tubes, ignoring corrections from more distant tubes and 
many-body effects associated with the screening of van 
der Waals interactions by the adjacent tubes Al- 
though these approximations introduce some errors in the 
values of E t they do not affect the key predictions of this 
work, which are sensitive to the variation of £^(R). The 
potential from each neighboring tube was derived with 
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the method of Stan et al (l9|. Because V(r,R) varies 
rapidly with r, values of E t (R) include large anharmonic 
and small anisotropic contributions. 

The form of N(E) for a given collection of nanotube 
bundles depends on sample preparation. A sample is rep- 
resented by an ensemble of points (one for each IC) in R 
space. The density of points in R space, a function /(R), 
is defined so that /(R) <iR is the number of IC's within 
an infinitesimal volume dR — dR\ di?2cW?3, centered at 
R. The function / enters the transverse density of states 
through this expression: 



g{E) = J dH f(R)6[E — E t (R) 



(2) 



At very low T, we need the value of /(R) only in 
the immediate vicinity of the global minimum of the 
energy (E rn at R TO ), but at higher T the specific ex- 
perimental distribution affects the results quantitatively. 
Here, we have assumed that the IC's are uniformly dis- 
tributed in R space within a radius spread of width 
3 A near R m . More general results will be described 
in a complete report of this work . We find, as one 
might expect, that R m occurs along the diagonal, sym- 
metry line (Ri = R2 = -R3). Along this line, which 
we call the (1,1,1) line, there occurs a global minimum 
E m = -1052.97 K at i?,=9.95 A. To derive g(E), it is 
important to know the variation of E t (R) near this mini- 
mum. Consider a particular displacement (in the Ri — R3 
plane) from R TO to a neighboring point for which R± has 
the same value, while R2 and R3 are slightly different: 
i?2 = Ri + S and i?3 = i?i — S. This change, parallel to 
the (0,1,-1) direction, yields an extremely small increase 
in the H2 energy, indicative of a very slow variation of the 
function E t (R) near R m . This behavior, shown in Fig.l 
reveals a long valley of low energy states in this direc- 
tion; the same behavior occurs along the five equivalent 
directions, e.g. (1,0,-1). 

The transverse density of states g(E), from Eq. 2, 
is shown in Fig. 2. Note that g{E) is proportional 
to \JE — E m for small E — E m ; the prefactor is deter- 
mined by the principal axes of curvature of the function 
E t (R). This square root behavior is identical to that 
found near a van Hove singularity in the phonon density 
of states of a 3D system near a minimum in the Brillouin 
zone, for the same reason- phase space topology [2(J ■ At 
higher energy, instead, the behavior of g(E) switches to 
1/V-E — E m . This (-1/2) power law follows from the fact 
that the higher energy displacements from R m are quasi- 
1D. As seen in Fig.l, the iso-energy contours are perpen- 
dicular to the diagonal, so that the energy gradient in 
R space is along the diagonal, with essentially constant 
transverse variation, a ID situation. 

Fig. 2 also presents N(E), derived with Eq. 1, which 
convolutes the transverse spectrum with the ID motion 
along the axis. The resulting power law behavior can be 
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FIG. 1: Contour plot showing the variation -Et(R) — E m near 
the minimum R m (which is chosen as the origin). From the 
contour closest to the minimum (the closed one) to the more 
distant ones, the energy contours correspond to Et(R) — E m = 
0.005, 0.01, 0.05, 0.1, 0.5, 1 and 5 K. 
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FIG. 2: Upper panel: Calculated transverse density of states 
g(E). Lower panel: Total density of states N(E), obtained 
from Eq.l in the text. E m = —1052.97 K is the lowest energy 
level corresponding to Ri — R2 = R3 — 9.95 A. 
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understood from realizing that if g(E) is proportional to 
(E — E m ) n , for some n, then N(E) is proportional to 
(E - E m ) n+1 / 2 . Hence, we find that N(E) is propor- 
tional to (E — E m ) near threshold. This linear behavior 
is that characteristic of a 4D gas in free space; the re- 
sult implies that this system exhibits J^D gas behavior at 
low T. For E > 20 mK above threshold, instead, N(E) 
becomes approximately constant, corresponding to the 
density of states of a 2D gas. Thus, there arises a di- 
mensionality crossover originating from the anomalous 
transverse density of states. We emphasize that the 4D 
regime is a direct consequence of the existence of a mini- 
mum in the function E t (R) , a result that is not sensitive 
to the details of the calculation [2l| ■ 

The thermodynamic behavior of the system is derived 
with the usual bose gas theory. For a given number 
of molecules, N, the chemical potential /i is determined 
from the relation 
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FIG. 3: Density of Hh molecules as a function of BEC transi- 
tion temperature. 



N=[ d E *W (3) 

At a specified T = l/(fcg/3), this relation yields a max- 
imum value N max when /i equals the lowest energy of 
the system, E m , at which point BEC begins. That is, a 
macroscopic fraction of the molecules fall into the low- 
est energy state when N > Nmax- Equivalently, at fixed 
N, BEC occurs when T falls below the inverse function 
T c = T{N max ). The resulting dependence on N of T c 
is shown in Fig. 3. As seen there, T c is of order 10 
mK, which is experimentally accessible. Fig. 4 shows 
the specific heat Cjv(T), calculated from the energy as 
a function of [i and T. The novel behavior observed in 
the figure is a result of the unusual form of N(E). At 
low T, Cn is proportional to T 2 , a consequence of the 
4D (linear) variation of N(E) at low E. Note the pres- 
ence of singular behavior of Cn(T) as one approaches 
the transition from above and a cusp at T c itself. At 
relatively high T, Cn(T)/ (Nks) is essentially unity be- 
cause N(E) is 2D-like at high E; a nondegenerate 2D 
gas has C N {T)/(Nk B ) = 1. Behavior for T > 0.2 K, not 
shown, is very sensitive to the distribution of nanotube 
radii. 

Discussions of these results with colleagues have led to 
several questions concerning the transition. One is this: 
since each IC contributes a density of states characteris- 
tic of a ID system, why does non-ID behavior arise here? 
The answer is that particles can exchange between dif- 
ferent IC's because of their common access to the vapor 
region. In practice, this may be a very slow process, lead- 
ing to nonequilibrium behavior. The presence of breaks 
or holes in the tubes may alleviate this kinetic problem 
without invalidating the model, since heterogeneity is 
an essential aspect of the model. Another question is 
whether similar behavior occurs for other bose gases. In- 
deed, BEC of 4 He is predicted by a similar analysis to 
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FIG. 4: Temperature dependence of the specific heat (per 
molecule) for molecular density n = 1 X 10~ 5 A -3 . 

occur if the sample's distribution of tubes includes those 
with some near its energy minima in R space, which 
occur near 8.5 A. Finally, one might wonder about the 
effects of intcrparticle interactions, which have been ig- 
nored up to this point |22| . Indeed, some previous studies 
of H2 in IC's have found that a liquid- vapor condensation 
occurs at a higher temperature (ss K) (in the absence of 
heterogeneity) 11]. A very recent study, however, found 
that nanotubes' screening of the intermolecular interac- 
tion reduces T c to about 10 mK |l8|. However, that 
calculation omitted the role of heterogeneity, which is 
relevant, according to the Harris criterion, since the 3D 
specific heat critical exponent is positive |23j . We expect 
that this condensation temperature is further reduced by 
disorder, enabling the BEC transition to occur. 

We summarize our results as follows. Heterogeneity al- 
ters the qualitative behavior of the low energy spectrum 
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of H2 molecules. The lowest-lying states of the system 
are those of particles in that channel. As T falls, particles 
aggregate in the (essentially 4D) space of quantum states, 
(R,p), with bose statistics having a dramatic effect, i.e. 
BEC, below a transition temperature of order 20 mK. 
Anomalous behavior is predicted for the specific heat, a 
consequence of the unusual density of states, which is 4D- 
like at very low energy and 2D-like at somewhat higher 
energy. An experimental probe of the real-space molec- 
ular density should reveal the needle-like concentration, 
below T c , of a macroscopic fraction of the particles within 
the lowest energy channel. 

Most intriguing to us is that this transition is a di- 
rect consequence of disorder, since the perfectly uniform 
system of identical nanotubes yields strictly ID, nonsin- 
gular behavior. Such a dramatic effect of heterogeneity 
occurs elsewhere in low temperature physics. Examples 
include the spin-glass transition [24j], the effect of tun- 
neling states on thermal behavior of glasses [2^] and the 
effects of disorder on monolayer films |2r3 ]. 
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